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Проведено численное исследование сопряжённого теплообмена и гидродинамики при движении вязкой жидкости в открытой 
прямоугольной полости с учётом процесса её разрушения. Получены гидродинамические картины вязкого течения в полости в 
условиях её разрушения и температурные поля для твёрдой и жидкой фаз. Изучено влияние динамических параметров и гео- 
метрических характеристик на формирование контуров образующихся кратеров. 


В последние годы наблюдается повышенный 
интерес исследователей к изучению конвективных 
течений в полостях различных типов. Большое ко- 
личество всевозможных теплотехнических уст- 
ройств, используемых в энергетических установ- 
ках, производственных процессах разного уровня 
сложности таких отраслей промышленности, как 
энергетическая, горная, космическая и многих 
других, представляют собой контуры открытого 
или закрытого типов, перенос тепла в которых осу- 
ществляется конвективной циркуляцией заполня- 
ющей их жидкости [1-5]. 

Процессы, связанные с изменением фазового 
состояния при плавлении материала протекают в 
условиях аварийных ситуаций, возникающих на 
энергетических установках по мере развития опас- 
нейших производств и усложнения технологичес- 
ких циклов в таких отраслях промышленности, как 
атомная, химическая, металлургическая [6-10]. 
Разработка оптимальных технологических режимов 
и выбор соответствующих параметров, обеспечива- 
ющих безопасность установок в условиях протека- 
ния производственных процессов, сопряжена с 
большими затратами времени и средств из-за реа- 
лизации в основном эмпирического подхода к ре- 
шению возникающих проблем данного типа. Кро- 
ме того, опытная отработка технологий всегда пред- 
полагает высокую вероятность выбора не самого 
оптимального режима, который при дополнитель- 
ных термических нагрузках, вызванных чрезвычай- 
ными ситуациями на объектах, может привести к 
авариям [6, 8, 9]. Возникает необходимость в теоре- 
тическом исследовании основных закономернос- 
тей гидродинамики и тепломассопереноса при дви- 
жении вязкой несжимаемой жидкости в полостях с 
изменяющейся геометрией области течения. 

Рассматривается нестационарное взаимодей- 
ствие ламинарной вязкой струи несжимаемой неи- 
зотермической жидкости с открытой прямоуголь- 
ной полостью [11]. Цель данной работы - исследо- 
вать процесс движения вязкой несжимаемой неи- 
зотермической жидкости в открытой полости пря- 
моугольного типа в условиях сопряжённого тепло- 
обмена с учётом изменения фазового состояния 
материала выемки (плавления). 


Исследования описанных процессов проводи- 
лись с использованием математической модели на 
основе системы уравнений Навье-Стокса в пере- 
менных вихрь-функция тока при умеренных чис- 
лах Рейнольдса 100<Ке<800, уравнений энергии 
для материала прямоугольной полости и жидкой 
фазы с соответствующими краевыми условиями, 
представленными в работе [11]. 

При исследовании движения вязкой жидкости 
в прямоугольной полости с учётом процесса плав- 
ления определяется нестационарное двумерное 
температурное поле, с использованием условия 

Ѳ > Ѳ, (1) 


определяется фронт плавления материала выемки. 
В формуле (1) 0* - безразмерная температура плав- 
ления материала. Условия на границе раздела твёр- 
дой и жидкой фаз в отличие от постановки [11] 
примут следующий вид: 

- на границе раздела между фазами в направле- 
нии по оси х выставляется граничное условие 
четвёртого рода для уравнения энергии, где вво- 
дится тепловой эффект плавления материала 
полости 
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- аналогично на границе раздела между фазами в 
направлении по оси у выставляется граничное 
условие такого же типа 
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В формулах (2, 3) соответственно приняты сле- 
дующие обозначения: 0 - тепловой эффект плав- 
ления материала, Я ь Я — коэффициенты теплопро- 
водности для твёрдой и жидкой фаз соответствен- 
но, р - плотность, ѵ - кинематическая вязкость, 
Г* - температура плавления материала полости, 
Т\ - начальная температура стенок полости, 
Рг - число Прандтля. 
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При этом скорости движения границы раздела 
между твёрдой и жидкой фазами в условиях изме- 
нения фазового состояния материала полости 
вдоль осей х, у соответственно будут определяться 
следующим образом 

у х і (г + Аі)-х,(г) ' у у,- (Т + А?) -у,- (?) 

1 Аі ’ 2 А( 

где х іг уі - координаты узлов разностной сетки на 
границе раздела фаз по осям пространственных пе- 
ременных х и у в моменты времени (+А( и і соотве- 
тственно, Ѵі, Ѵ 2 - скорости движения границы раз- 
дела фаз в условиях изменения фазового состояния 
материала полости в направлении по осям прост- 
ранственных переменных хи у. 

При исследовании данного процесса необходи- 
мо выполнение интегрального соотношения, явля- 
ющегося основной интегральной характеристикой 
течения: 

х 0 х 1 
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здесь X] - координата фиксированной крайней точ- 
ки входного участка, лежащей на линии симметрии; 
х« - координата точки раздела между участками с 
разным направлением движения жидкости в полос- 
ти; х 2 - координата крайней точки кратера, образую- 
щегося при плавлении материала полости; х 3 - ко- 
ордината фиксированной крайней точки выходного 
участка, лежащей на боковой стенке полости; 
ѵ_(х,у) - поперечная составляющая скорости дви- 
жения жидкости в направлении от входного участка 
до дна полости; ѵ_|_(х,у) - поперечная составляющая 
скорости движения жидкости в направлении от дна 
полости до выходного участка; Щ - расход жидкос- 
ти за счёт плавления материала выемки. 

Постановка задачи сделана при следующих до- 
пущениях: 1) двумерность; 2) материал жидкой и 
твёрдой фаз идентичен; 2) теплофизические харак- 
теристики материалов твёрдой и жидкой фаз при- 
няты идентичными; 3) реализуется ламинарный 
режим течения; 4) возможные эффекты газовыде- 
ления, обусловленные химическими и физически- 
ми процессами, не учитываются. 

Методом конечных разностей решены система 
уравнений Навье-Стокса в переменных функция 
тока-вихрь, уравнение энергии и уравнение теп- 
лопроводности с соответствующими граничными 
условиями. Разностные аналоги уравнений пере- 
носа и теплопроводности решены методом прогон- 
ки [12-14]. Уравнение Пуассона на каждом вре- 
менном слое решалось методом последовательной 
верхней релаксации. Граница фронта плавления 
определялась с помощью критериального условия 
фазовых переходов (1). При несовпадении криво- 


линейной границы в условиях изменения фазового 
состояния материала полости с узлами сетки при- 
менялись разности с переменными шагами по ко- 
ординатам. Использовалась разностная схема вто- 
рого порядка точности. Были использованы равно- 
мерные и неравномерные разностные сетки. Зна- 
чения вихря на границе гидродинамической облас- 
ти, а также на подвижной её части, определялись из 
уравнения для функции тока. 

Изучение гидродинамики и сопряжённого теп- 
лообмена при движения вязкой жидкости в полос- 
ти прямоугольного типа было проведено в услови- 
ях изменения фазового состояния материала выем- 
ки (плавления). Фронт плавления материала по- 
лости определялся с помощью полученного дву- 
мерного неоднородного поля температур с исполь- 
зованием критерия (1). Течение в условиях описан- 
ного процесса изучено в широком диапазоне изме- 
нения динамического параметра 100<Ке<800 и 
других параметров модели. В результате взаимо- 
действия высокотемпературной струи вязкой жид- 
кости с прямоугольной полостью наблюдается об- 
разование кратеров параболической формы, опре- 
деляющих фронт фазового перехода материала. Ти- 
пичные результаты численных исследований при- 
ведены на рис. 1-3. 

На рис. 1, а представлены контуры кратеров, 
образующиеся в результате плавления полости 
прямоугольного типа в различные моменты време- 
ни 1, 875-ЮЛ 2,75- ЮЛ 3,75-10"*. Согласно анализу 
приведённого рисунка можно отметить, что кри- 
вые 1, 2, 3, определяющие контуры кратеров, в 
симметричной области 0,45<х<0,50, 0,11 <>»<0,25 
более вытянуты вдоль линии симметрии по сравне- 
нию с полушириной образованной лунки. С увели- 
чением времени скорость плавления стенок полос- 
ти вдоль оси х существенно уменьшается по срав- 
нению со скоростью плавления вдоль оси у - осо- 
бенно отчётливо эта тенденция наблюдается на 
контурах 2, 3. На рис. 1-3 обозначено: 4 - линия 
симметрии; 5 - граница раздела фаз. 

На рис. 1, б представлены контуры кратеров, 
образующиеся в результате плавления полости 
прямоугольного типа в различные моменты време- 
ни 1,0-ЮЛ 1,875-10"* при а = 0,32а* (а - длина 
входного участка, а* - длина проницаемого участ- 
ка полости). При указанных параметрах размеры 
осесимметричной области, в которой происходит 
кратерообразование, увеличиваются 0,43<х<0,50, 
0,095<у<0,25. Следует отметить, что при данных 
параметрах скорость плавления полости вдоль оси 
у значительно превосходит скорость плавления 
вдоль ОСИ X. 

На рис. 1, в показаны контуры кратеров, обра- 
зующиеся в результате плавления полости прямоу- 
гольного типа в различные моменты времени 
1,5-ЮЛ 2,375-ЮЛ 3,375-10"* при а = 0,16а*. 

На рис. 1, г приведены контуры границ гидро- 
динамической области в условиях изменения фазо- 
вого состояния вещества при плавлении полости в 
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различные моменты времени 1,551 -10^, 2,426* 10^*, 
3,426- 10 -4 при а = 0,24а*. Сравнение рис. 1, 6-г от- 
ражает динамику изменения контуров кратеров в 
условиях плавления вещества при определённых 
выше параметрах. 

При движении жидкости от нижней границы 
гидродинамической области к выходному участку 
(второй этап движения жидкости) поперечная сос- 
тавляющая скорости начинает расти с увеличением 
координаты по оси пространственной переменной 
у, а полуширина сформировавшегося контура кра- 
тера постепенно возрастает, при этом максималь- 
ное значение её меньше высоты образующейся 


лунки. Достаточно хорошо это просматривается на 
графиках (рис. 1, а-г), отражающих формировние 
контура фронта плавления. За счёт нагрева матери- 
ала полости, обусловленного теплоотводом в её 
стенки, изотерма плавления сдвигается по отноше- 
нию к начальному её положению, в результате чего 
происходит искривление контура фронта фазового 
перехода материала выемки. 

Проведённый анализ позволяет отметить влия- 
ние входного параметра а при одних и тех же зна- 
чениях других параметров модели на процесс кра- 
терообразования при плавлении стенок полости. С 
увеличением динамического параметра и умень- 





Рис. 1 . Контуры кратеров, образующиеся в результате плавления полости в различные моменты времени при следующих па- 
раметрах: а) Ре=400, Рг=0,979 ; б) Ре=600, Рт=0,72 ; в) Ре=600, Рт=0,979; г) Ре=600, Рг=0,979 
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шением длины входного участка качественная сто- 
рона кратерообразования существенно не меняет- 
ся. Следует только отметить, что кратерообразова- 
ние начинается раньше, и уменьшаются размеры 
симметричной области, в которой протекает дан- 
ный процесс: 0,46<х<0,50, 0,11<у<0,25. 



Рис. 2. Распределения температуры в твёрдой и жидкой фа- 
зах в условиях плавления полости в различных сече- 
ниях у по координате х в фиксированный момент 
времени при числах Ре=400, Рт=0,979 

На рис. 2 приведены распределения температу- 
ры в твёрдой и жидкой фазах в условиях плавления 
в фиксированный момент времени в различных се- 
чениях у: 0,216, 0,256, 0,302. Как видно из рисунка, 
температурный профиль до границы раздела фаз 
достаточно быстро убывает. На границе раздела 
фаз возникают некоторые локальные минимумы 
температуры, что можно объяснить тепловым эф- 
фектом плавления. Изменения температурного 
профиля в твёрдой фазе носят менее выраженный 
характер, чем в жидкой (кривые 1, 2, 3). Это можно 
объяснить тем, что нагрев твёрдой фазы происхо- 
дит лишь за счёт теплопроводности, в то время как 
охлаждение жидкой фазы происходит в условиях 
конвективного переноса. 

На рис. 3 представлены профили поперечной 
составляющей скорости ѵ в фиксированный момент 
времени внутри кратера, образующегося в результа- 
те плавления стенок полости прямоугольного типа. 
Результаты расчёта приведены для сечений у: 0,096, 
0,376 при а =0 ,34а*. На первом этапе [11] движения 
вязкой несжимаемой неизотермической жидкости 
профиль поперечной составляющей скорости убы- 
вает по всей высоте образующегося кратера. На вто- 
ром этапе движения жидкости профиль поперечной 
составляющей скорости начинает расти по всей вы- 
соте кратера, при этом профиль возвратного тече- 
ния претерпевает деформацию за счёт сложного 
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контура кратера. Степень деформации профиля по- 
перечной составляющей скорости возрастает с 
уменьшением высоты кратера (кривые 1, 2). 



Рис. 3. Твёрдая фаза (заштрихована) и профили попереч- 
ной составляющей скорости ѵ в различных сечениях у 
при Ре=600, Рт=0,979 


На основании полученных результатов можно 
сделать вывод о том, что характер поведения тем- 
пературных профилей в различных сечениях по- 
лости и оригинальных контуров кратеров качест- 
венно сохраняется с изменением динамического 
параметра Яе, претерпевая незначительное изме- 
нение при варьировании параметра входа а. 

Результаты, полученные при исследовании 
описанного процесса, позволяют проследить дина- 
мику процесса плавления полости, показать влия- 
ние параметров модели на процесс кратерообразо- 
вания и условия теплообмена при движении жид- 
кости в открытой выемке. 

Анализ проведённых исследований позволяет 
сделать вывод о возможности дальнейшего расши- 
рения области применения сформулированной ма- 
тематической модели для решения задач о конвек- 
тивных течениях в полостях различных типов в ус- 
ловиях изменяющейся геометрии области течения 
жидкости. В дальнейшем можно решать задачи та- 
кого типа для областей со значительно более слож- 
ной геометрией (условия криволинейных контуров 
расчётных областей), чем рассмотренные в данной 
работе и в [14, 15]. 
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